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We investigate the isotropic-to-nematic phase transition in systems of hard helical 



"^3 particles, using Onsager theory and Monte Carlo computer simulations. Motivation 

O of this work resides in the ubiquity of the helical shape motif in many natural and syn- 



thetic polymers, as well as in the well known importance that the details of size and 



K*- shape have in determining the phase behaviour and properties of (soft) condensed 

^O matter systems. We discuss the differences with the corresponding spherocylinder 

^~~. phase diagram and find that the helix parameters affect the phase behaviour and 
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I. INTRODUCTION 

In the physics of fluids, simple and complex, it is well established that size and shape of 
particles, either molecular or colloidal, play a key role in determining thermodynamics, 
structure and dynamics (e.g. Refs. [HE])- Hard body particles, interacting with each 
other through steep repulsive potentials only, can thus be viewed as elementary models to 
understand the phase behaviour and properties of physical systems. 

While the first example of phase transition driven by purely steric interactions is un- 
doubtely the fluid-to-crystal phase transition in hard spheres Pj, the isotropic-to-nematic 
liquid crystal (IN) phase transition in hard slender rods predicted by Onsager |1] paved the 
way for an entirely new field. Although Onsager theory was originally motivated by the ob- 
servation of a nematic liquid crystal phase in suspensions of inorganic and biological rod-like 
colloidal particles |5j, its influence over the years have proven to be much more profund. The 
explanation of an ordered fluid formation, such as the nematic phase, as the result of the 
competition between orientational and translational entropy contributions, introduced the 
far reaching concept of "ordering entropy", whilst the idea of expressing the system's free 
energy as a functional of single particle density is a precursor of modern density functional 
theory (DFT) (e.g. Ref. [6]). 

The original Onsager theory accounted only for the second-virial coefficient contribution, 
thus (strictly) limiting its applicability to rod-like particle systems with large aspect ratios, 
but several improvements have been more recently proposed to overcome this drawback 
and include also higher order contributions. This prompted a number of DFT approaches 
with different degrees of sophistication, as well as a series of computer simulations, [3-[9], 
that can be applied and extended to many systems, either mono- or poly-disperse, both 
homo- and hetero-geneous. All these studies have confirmed that entropy alone can lead 
to complex phase organisations, including smectic and columnar liquid crystals. Recent 
simulations have further unveiled examples of the zoo of morphologies that can be obtained 
from packing of particles of different shape fiU\ [TT] . 

So far, most theoretical and computational studies have focussed on convex hard particles 
(e.g. Ref. [7j), whilst concave particles have been given less attention. Besides simple 
dumbbells, these include bent-core p^ [13] . lens-like [T^ and bowl-shaped [15] particles. 
Somewhat surprisingly, hard helices are not part of the above list, in spite of the several 



examples of this shape that can be found in natural and synthetic polymers. Rigid and 
semiflexible helical polymers (polynucleotides, polypeptides, viruses) have a well known 
propensity to form liquid crystal phases at high concentration [TBHT^ . When examining and 
interpreting the experimental phase behaviour, helicoidal particles were generally assimilated 
to rods, thus neglecting peculiarities related to the actual shape (e.g. [20]). 

In this work, we address this problem by undertaking a study of the phase behaviour of 
hard helices as a function of their structural features. The model helices, shown in Fig. [T| are 
obtained by considering a set of fused hard spheres all having diameter D, and arranged in a 
helical fashion along a string of contour length L. Helices of different shapes and Euclidean 
lengths A are generated by tuning the radius r and the pitch p (the helix parameters are 
defined in Appendix A). 
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FIG. 1: Model helix and its characteristic parameters. 



Considering systems of rigid homochiral helices, we have focussed on the effects of the 
particle shape on the IN phase transition. Homochiral helices are expected to form a chiral 
nematic phase (cholesteric), in which the local preferred orientational axis (the director) 
rotates in space around a perpendicular axis. However the cholesteric pitch is known to 
be orders of magnitude larger than length scale of inter-particle interactions; thus the local 
phase properties are virtually indistinguishable from those of the corresponding untwisted 
nematic phase and the influence of chirality on the phase boundaries can be neglected at this 
stage [21]. In addition, the presence of periodic boundary conditions in the MC computations 



does not allow the emergence of an equilibrium cliolesteric order, with a pitch much longer 
than the size of the simulation box. For these reasons, phase chirality will not be considered 
in the present study and will be the object of a future dedicated study. 

Our approach hinges upon an Onsager theory supported by Monte Carlo (MC) simula- 
tions. We have used a simple extension of the Onsager theory suitable for particles with a 
finite aspect ratio, where the role of higher order virial terms is effectively taken into account 
[22t |23] . This theory was previously applied with success to the thermodynamics of the IN 
phase transition in hard, straight rods [23] • Here its performance for the more subtle and 
challenging case of hard helices has been tested by comparison with MC simulations. As 
further elaborated below, MC simulations for hard helices are considerably more demand- 
ing with respect to the spherocylinder counterpart. Onsager theory, being computationally 
much less expensive, can then be rather useful for a preliminary exploration of the phase 
diagram. 

Through a combined action of theory and simulations, we will provide evidence that the 
helical shape significantly affects the location of the IN transition in comparison with hard 
rods of similar length and diameter, thus casting some doubts on the possibility of providing 
a straightforward link between the phase behaviour of hard helices and that of hard rods. 

The layout of this work is as follows. In the next section our Onsager and MC approaches 
are illustrated, while theoretical and simulation results are presented and discussed in Section 
III. Section IV concludes this work by summarising present findings and providing an outlook 
for future developments. 

II. THEORY AND SIMULATIONS 

A. Theory 

1. Free energy of the isotropic and nematic phases 

Let us consider a system of A^ identical helices in a volume V at temperature T. We denote 
hj V = V/N the volume per particle and rj = vq/v the packing density, where Vq is the 
volume of a particle. 



The mutual interaction between a pair of hard hehces (1 and 2) takes the form: 

{oo if 1,2 overlap 
(1) 
if 1,2 do not overlap 

where -Ri and R2 are the positions of the center-of-mass for helices 1 and 2 respectively, 
-R12 = -R2 — -Ri is a vector defining the relative position of helix 2 with respect to helix 1 
and f2i2 = (ai25 /3i2,7i2) are the Euler angles that define the rotation from 1 to 2. 

In the Onsager approach the free energy of the system is expressed as a functional of the 
single particle density function p{R, ft), where R is the particle position and ft = (a, (3, 7) 
are the Euler angles specifying the particle orientation, with the normalization condition 
J dRdflp{R, fl) = N. The nematic phase is uniform and the density function depends only 
on the particle orientation. Assuming helices as uniaxal particles, the single particle density 
in the uniaxial nematic phase (see below) reduces to |24J: p = p{f3) = /(/3)/(47r^t>), where 
P is the angle between the helix axis and the nematic director and /(/?) the orientational 
distribution function. The latter obeys the normalization condition J_^d{cos /3)f{f3) = 1 
and is a constant equal to 1/2 in the isotropic phase. 

The free energy can then be expressed as: 



A[f{/3)] = NksT 



A^ B 
In^^ + lnAT-l 



+ A-[/(/3)] + A-[/(/3)] (2) 



The first term is the Helmholtz free energy of the ideal gas: Atj. = {h? /2'KkBTm) is the 
de Broglie wavelength and Gor = h'^/Sn'^kBl is the rotational temperature, with fc^ and h 
being the Boltzmann and the Planck constant, respectively, while m is the mass and I is 
the inertia moment of the particle. The second term in eq. ^ accounts for the decrease of 
orientational entropy due to orientational ordering: 



NksT J ^d{cosP)f{P)ln[2fm (3) 

and the last term, A^^, represents the excess free energy. Within the Onsager formulation 
coupled with the Parsons-Lee (PL) correction [221 123], this is expressed in terms of the 
second virial contribution with a pre-factor, G{r]), that is meant to account for higher virial 

contributions: 

^^^ G{ri) 



NksT 2(47r2)2i; 



J dn^ /(/3i) J dn^ f{/32)ve.ci{nu) (4) 



where ^6x01(^12) is the volume excluded to a helix (2) by another (1): 

'^6x01(^^12) = - / dRi2ei2{Ri2,^i2), (5) 

with the Mayer function [2j: 

euiRi2, f^i2) = exp{-U{Ru, nu)/kBT} - 1. (6) 

Introducing the second virial coefficient: 

B2 = ^j^ j dn^f{f3,) j dfl2f{(i2)VeU^12). (7) 

the excess free energy eq. Q becomes 

A^^[f]/NkBT = G{v)B2[f]/v. (8) 

2. Parsons-Lee (PL) and Modified Parsons-Lee (MPL) approximations 

The approximation proposed by Parsons [22] and subsequently used by Lee |2S! and others 
[23 [2S] for hard sperocylinders (and ellipsoids) relies on the assumption that the excess free 
energy is proportional to that of a system of hard spheres (HS) at the same packing fraction 

^-(r^) _ A-5(r/) ^g^ 



NksTB^iv) NkBTB^^iji) 
Use of the Carnahan-Starling expression for the free energy of hard spheres [2Z], along with 
the relationships B^^ = Avhs and rj = vhs/v, where vhs is the volume of a hard sphere, 
yields 

GM = -^^^^^^^ = 11^^ (10) 

In the original and subsequent works [23| |25| 126] the volume of the reference hard spheres 
was taken equal to that of the spherocylinders (or ellipsoids), vhs = "^o- Good agreement 
between theory and simulations was obtained in that case, but significant discrepancies were 
found for linear particles made of tangentially bonded hard spheres |28j. It has then been 
suggested that the assumption vhs = "^o may be inappropriate for hard non-convex bodies, 
since in this case the free volume available at a given number density is smaller than for 
convex particles having the same geometrical volume [2^ . It was proposed that in this case 
the volume of the reference hard spheres should be replaced by an effective volume, Vei, 



defined as the volume of tlie non-convex particle that is inaccessible to other particles. This 
effective volume is larger than the geometrical volume, and for linear hard spheres it has 
been evaluated in ref. [30] (see Appendix IbI). This variant of the PL theory has been given 
the name of modified Parsons-Lee (MPL) theory |29j . 

3. Expansion in terms of orientational order parameters 

The orientational distribution function /(/3) is conveniently expanded in a basis of Legendre 
polynomials 

^ oo 

/(/5) = 2 E (4J + 1) (^2,) P2,(cos/3) (11) 

j=0 

where {P2j) are the nematic order parameters 

(P2,) = j rf(cos/3)/(/3)P2,(cos/3) (12) 

The expansion is limited to polynomials of even rank in view of the nonpolar character of 
the nematic phase. The order parameters take values in the range —1/2 < {P2j) < 1 and 
vanish in the isotropic phase. 

Upon substituting Eq. ( 11 ) in Eqs. ([sJ-Q and exploiting the properties of Wigner rotation 



matrices [3l], we can express the orientational and excess contributions to the Helmholtz 
free energy as a function of the order parameters: 



NkeT 2 . ^ 



1 oo „i 

- J2 (4J + 1) (^2,) / rf(cos /3)P2,(cos /3) In 



^(4/ + l)(P2,,)^2,'(cos/3) 



j'=0 



(13) 



A*'^ G{r]) 



oo „ 

J](4j + 1) {P2,f I dn^2 P2j{cOS/3i2)Ve.cl{nu) (14) 



" 3=0 

This leads to the following expression for the pressure 

P 1 fdA\ 1 1 /^, s r/(5-3r/) 



keT kBT\dVj^j. v lOvrV V 4(1 - r/)3 

oo „ 

5^(4j + 1)(P2,)' / dn,2P2AcOS/3,2)Ve.cl{^l2) 



OO 

X 

j=0 



(15) 



B. Monte Carlo simulations 

In order to test the theoretical predictions, we implemented Isothermal-Isobaric (NPT) 
MC simulations j32l 133] on a system of A^ hard helices, contained in cubic or orthorhom- 
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bic computational boxes, with the usual periodic boundary conditions. Simulations were 
organised in cycles, each consisting of 2N attempted particle moves (a random translation 
and rotation) and a volume move. Rotation trial moves were implemented either using the 
Barker- Watts [M] or the quaternions methods [351 ES], finding a good consistency between 
them. Volume moves were either performed scaling up or down the box in those cases where 
cubic boxes were used or attempting to change a randomly selected edge of the box in the 
other cases. Being concerned with the IN phase transition only, we neglected from the outset 
all possible nuances necessary to properly account for other phases, such as twisted-nematic 
boundary conditions and variable-shape computational boxes. 

The overlap condition was computed by first inserting each helix into the smallest sphero- 
cylinder containing it and testing for overlap between two such spherocylinders (see FigM. 
This is a relatively fast test as it amounts to finding the minimal distance between two 
segments. 

To this purpose, we used the algorithm proposed by Vega and Lago [2Z]. This method is 
approximatively four times faster than others previously used, essentially because it reduces 
to only four the number of regions to be checked for closest approach. 

Only in the event of overlap between two spherocylinders, the spheres forming the em- 
bedded helices were tested for overlapping. This procedure significantly reduced the com- 
putational cost of the overlap test, that is one of the bottlenecks of this type of simulations, 
and considerably increased their efficiency. 




FIG. 2: Cartoon of the overlap between the spherocyhnders containg a pair of hehces. 
The IN phase transition was monitored using the main orientational order parameter. 



(P2), already discussed in the theoretical section. To this aim, the following tensor [38] 

1 ^ S 1 

was evaluated, and the corresponding eigenvalues and eigenvectors computed. Here a, /3 = 
X, y, z and u" is the a component of the unit vector Uj describing the orientation of the 
2— th helix axis. The orientational order parameter (P2) was then identified with the largest 
eigenvalue of Q. The difference between the other eigenvalues of Q was found to be smaller 
than 5%, in agreement with our assuption of uniaxial nematic order. 

C. Computational details 

1. Onsager theory calculations 

For each system, the Helmholtz free energy given in eqs. ([2]), (13) and (14), is minimised at 
increasing values of the density 1/v, and the order parameters of the stable phase at each 



density value obtained are then used to calculate the pressure according to eq. (15). 



a. Evaluation of pair integrals Integrals over all the relative positions and orientations 



of pairs of particles, appearing in eqs. (14) and (15), are preliminarly evaluated and stored, 
to be used for the calculations at the various density values. These integrals have the general 
form: 



/ dai2 / rf(cos/3i2)P2,(/3i2) / ^712 / #12 / rf(cos^i2)(P?2)' (17) 

Jo J-1 Jo Jo J-1 

where R12, the vector position of molecule 2 with respect to molecule 1, is expressed in 

spherical coordinates, R12 = {-R12, 012,^12} and -R°2 is the closest approach distance, which 

is a function of the relative position and orientation of the two molecules. The computational 

cost of this sixfold integral scales with M^, where M is the number of spheres in a helix. 

Gauss-Legendre and Gauss- Chebyshev quadrature algorithms are used to evaluate these 

integrals |39]. 

b. Free energy minimization It is expedient to choose as variational parameters the 

coefficients Uj of the expansion 

00 
-ln/(/3) = ^M,P2,(cos/3), (18) 
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rather than the order parameters. Eq. (18) is used in eq. (tsl) for the orientational contri- 



bution to the free energy, A°^, and is introduced into eq. (14) for the excess contribution 



j^ex^ through the the order parameters, eq. (12). Thus, the Helmholtz free energy is ex- 



pressed as a function of the Uj coefficients. This has a twofold advantage: The expansion 



eq. (18) converges faster than that of the density function eq. (11) and the parameters Uj 



are unconstrained, unhke order parameters. 

2. Monte Carlo simulations 

Our NPT MC simulations were carried out using A^ = 675 or 867 hard helices with periodic 
boundary conditions. As a general rule, we started a series of simulations from a diluted 
configuration and reached equilibrium upon compression. Typical equilibration runs con- 
sisted of 3 X 10^ MC cycles and were followed by a production run of additional 3 x 10^ MC 
cycles, during which averages of density and order parameter were calculated. 

In most of the equilibration runs the maximum values for the displacement, rotation and 
volume moves were varied in the course of the run to reach a 30-40% of acceptance. This 
procedure is known to lead to a possible violation of the detailed balance condition |10], 
but we explicitly verified that this does not lead to any bias in the present case. During 
production runs, the overall acceptance ratio was adjusted to be 30% — 40% by a suitable 
choice of the maximum displacement, rotation and volume parameters, and these values 
were never altered during the run. 

It is worth emphasizing that simulations for hard helices are considerably more demanding 
from the computational point of view than simulations of hard spherocylinders. Depending 
on the state point considered and the values used for the radius and pitch, the computational 
cost might be as high as 8 times that of the corresponding spherocylinders. 

III. RESULTS AND DISCUSSION 

We have considered helical particles with different structural parameters (see Appendix lA] 
for the definition), as in Fig|3] Helices are formed by 15 fused spheres and have the same 
contour length L = lOD, but with different pitches and radii. 

In presenting and discussing our results we will use reduced units, with the diameter D 
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taken as the unit of length, and with reduced pressure P* = PD^/kBT. For each system, the 
results from MC simulations were compared with those from Onsager theory with the PL 
and the MPL approximation, which differ in the definition of the packing fraction entering 



the scaling factor G, eq. (10): rj = vq/v, with vq being the geometric volume of the helix 
(PL), and r] = Vef/v, where Ve/ is the effective volume defined in Appendix [B] (MPL). Values 
of geometric and effective volume are reported in Tables [T] and [TT| respectively. MC data 
will be reported with error bars, evaluated according to the reblocking algorithm described 
in ref. [H]. 

r=0.2 




p=2 p=4 p=8 
A=8.47 A=9.54 A=9.88 

r=0.4 




LHSC 

A=10 



p=2 p=4 p=8 
A=6.23 A=8.47 A=9.55 

FIG. 3: Helices of radius r and pitch p investigated here. Helices have the same contour length 
L but different Eucliedean lengths A (lengths are scaled with the the sphere diameter D). For 
comparison also the fully extended linear hard sphere chain (LHSC) is shown. 

As a preliminary test, we have performed calculations for the the linear hard sphere 
chain (LHSC), for which A=L. Figure 111 shows order parameter (P2) and reduced pressure 
P* calculated for the LHSC as a function of the packing fraction rj = v/vq. At 77 ~ 0.24 an IN 
phase transition occurs, characterized by a jump in the order parameter. On moving deeper 
in to the N phase (P2) takes higher values, larger than 0.8. The nonvanishing (P2) obtained in 
the isotropic phase from simulations can be attributed to finite-size effects, and this feature is 
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also present in the isotropic phase for hehces. Fig. |4] shows good agreement between theory 
and simulations for LHSCs. The results obtained using the PL and the MPL approximation 
are also very close one another, as expected in view of the high superimpositions of the 
spheres, so that the cavities between them have tiny volumes. This agrees with ref. [29] 
where it was shown that for LHSCs the discrepancies between MC simulations and PL 
theory, and correspondingly also the improvements of the MPL scaling, decrease as the 
superposition between adjacent spheres increases. 
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FIG. 4: {P'l) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction r] = vq/v for the LHSC, from MC simulations (closed circles) and from Onsager theory 
with PL (dashed line) or MPL (solid line) approximation. Insets on the left panel, here and in 
following figures, depict representative snapshots obtained using QMGA software. |42j 



Figures [5]j7] show order parameters and pressures calculated for the helices with r=0.2 and 
decreasing pitches p=8, 4 and 2. In all these cases a IN transition is clearly visible, with its 
location in densities shifting from rj ~ 0.24 to ?7 ~ 0.29 with decreasing pitches from LHSC 
(infinite pitch) to the helix with shorter pitch (p = 2). This can be qualitatively understood 
in terms of the decrease of the Euclidean length (and hence the aspect ratio) with decreasing 
pitch. In all these cases, we find a good agreement between Onsager theory and numerical 
simulations in the location of the IN transition and in the density dependence of the (P2) 
order parameter. However pressure tends to be understimated by theory, especially in the N 
phase, and this differences increase with increasing density and with decreasing pitch. The 
PL approximation does not appear to be adequate for these helical particles and use of the 
MPL variant leads only to a very slight improvement. The reason is that the non-convexity 



13 



of the helices is not simply due to the voids between adjacent spheres (see Fig. 12), so 
removal of these voids is not sufficient to account for the the real excluded volume. 
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FIG. 5: (^2) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction t] = vq/v for the helix with p= 8 and r= 0.2, from MC simulations (closed circles) and 
from Onsager theory with PL (dashed line) or MPL (solid line) approximation. 
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FIG. 6: (P2) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction t] = vq/v for the helix with p= 4 and r= 0.2, from MC simulations (closed circles) and 
from Onsager theory with PL (dashed line) or MPL (solid line) approximation. 



The discrepancy with respect to LSHCs becomes even more pronunced for larger radii, 
as depicted in Figures |8 10 reporting the (P2) order parameter and the reduced pressure 
calculated for helices with r=0.4 and p=8, 4 and 2. These helices are curlier than those with 
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FIG. 7: (P2) order parameter (Left) and reduced pressure P* (Right) as a function of the packing 
fraction r] = vq/v for the hehx with p=2 and r= 0.2, from MC simulations (closed circles) and 
from Onsager theory with PL (dashed line) or MPL (solid line) approximation. 



smaller r value (see Fig. pi), so it is not surprising that the differences from the behaviour 
of LSHCs are even more pronounced. No clear N phase is observed in simulations for the 
helices with p=4 and p=2, although at sufficiently high packing fraction (77 k. 0.35) an 
anisotropic organization, with some signature of layered ordering, is visible. A complete 
characterization of these phases is delicate, mainly due to equilibration problems, and is 
presently under scrutiny. In the case of p = 8 (see Figjs]), a nematic phase was detected 
between r] ~ 0.27 and t] ~ 0.38; interestingly, the IN transition occurs at higher density 
than for the helices with smaller radius and similar Euclidean length (r=0.2 and p=4). In 
all helices with r=0.4 we have also found a marked deviation between theoretical and MC 
results. In contrast to simulations, a nematic phase is predicted by Onsager theory for all 
pitch values, with the IN transition occurring at increasing density as the pitch decreases. Of 
course, being the theory implemented only for isotropic and uniaxial nematic phases, other 
possible phases could not be investigated. In short, only for the most elongated system 
(r=0.4 and p=8), we find a reasonable agreement between theory and simulations in this 
case. For shorter pitches, a jump in the (P2) order parameter is obtained from simulations 
and theory at similar 77 values, but the ordered phases appear to be different. As for pressure, 
differences between theory and simulations even appear in the isotropic phase for the helices 
with p=A and p=2, with theoretical predictions lower than the MC results. For p=2 the 
improvement deriving from the MPL approximation is more significant than in the other 
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cases, due to the larger value of the effective volume determined for this system using the 
rolling sphere criterion (see Table III] and Appendix IbI) . 





FIG. 8: (P2) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction ry = vq/v for the heUx with p= 8 and r= 0.4, from MC simulations (closed circles) and 
from Onsager theory with PL (dashed line) or MPL (solid line) approximation. 
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FIG. 9: (P2) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction r/ = vq/v for the helix with p= 4 and r= 0.4, from MC simulations (circles) and from 
Onsager theory with PL (dashed line) or MPL (solid line) approximation. Open circles are used 
for metastable states not yet fully characterized. 



An interesting last point, related to the above findings, is whether the IN phase transition 
for helices can be mapped on to that of rods in terms of simple parameters like the aspect 
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FIG. 10: (P2) order parameter (Left) and reduced pressure P* (Right) as a function of the volume 
fraction r/ = v^jv for the hehx with p= 2 and r= 0.4, from MC simulations (circles) and from 
Onsager theory with PL (dashed line) or MPL (solid line) approximation. Open circles are used 
for metastable states not yet fully characterized. 



ratio, as generally done in experimental work on helical systems |20]. Figure 11 collects 
the theoretical predictions of the IN phase transition as a function the Euclidean length 
A. For comparison, the results obtained for LHSCs and those for spherocylinders are also 
reported. In the latter case the Onsager expression for the excluded volume was used [1]. 
Of course the contour length L, which is identical for all the helices, is not a significant 



parameter in relation to the IN phase transition. On the other hand. Figure 11 also suggests 
that the Euclidean length, although more meaningful, is not fully satisfactory either, since 
for the same aspect ratio A/D, the density at which the IN transition occurs has a non 
trivial dependence on the combination of the helical parameters r and p. As a general rule, 
we find the transition to move towards higher volume fraction with increasing degree of 
non-convexity. The fact that the location of the IN phase transition is not uniquely related 
to the aspect ratio may have implications for the analysis of experimental data for helical 
particles, as anticipated. 



IV. CONCLUSIONS 



In this work, a system of hard helical particles have been investigated using an Onsager 
theory and MC computer simulations. Our main goal was to rationalize the changes in 
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the isotropic-to-nematic phase transition on going from straight rod-like to quite tortuous 
helical particles. We have found that helicity affects the location of the IN phase transition, 
with the latter in general being shifted to higher densities with increasing aspect ratios, as 
in spherocylinders. However, the same aspect ratio can be achieved with different struc- 
tural parameters of hard helices, and this affects the IN phase transition. In other words, 
the aspect ratio alone cannot be considered as a good candidate for the interpretation of 
liquid-crystal phase diagrams of strongly curled helical particles. Our numerical results also 
unveiled the presence of additional ordered phases, especially in the case of highly distorted 
partices, that will require further analysis to assess their specific natures. 

Another objective of our study was a test of Onsager theory for helical systems. We have 
examined whether the Onsager theory, which has proved to be successful in accounting for 
the thermodynamics of the IN phase transition in hard rods, can be extended to systems 
of helical non-convex particles. We have found that for high helicity Onsager theory signif- 
icantly departs from numerical simulations, even when a modifed form of the Parsons-Lee 
rescaling is included to account for the non-convexity of particles. When compared to the 
MC simulation data, Onsager theory generally understimates pressure, with deviations that 
increase with increasing density and upon going from the isotropic to the nematic phase. 
This points to the need of a more effective theory for hard non-convex particles, a field 
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that remains largely unexplored. Besides the Onsager theory employed here, various other 
theoretical approaches have been proposed, which include scaled-particle theory [iil - HB] , the 
Vega and Lago theory that aims at incorporating a better description of the isotropic state 
|47j . as well as Wertheim statistical mechanical treatment of associating fluids [48j, which 
was successfully applied to bent-shaped particles in the isotropic phase |26J. Another ap- 
proach envisages the extension of the Onsager theory beyond the second virial contribution. 
We intend to undertake a thorough analysis and comparison of available theories in a future 
work. 

This study is preliminary from several viewpoints and we plan to extend it in a number 
of ways, as alluded before. Though we dealt with chiral particles, phase chirality was not 
considered here, and the relationship between particle and phase chirality is one of the 
high priority points in our agenda. Other points that deserve close attention, and are 
currently under investigation, are a detailed numerical definition of the IN coexistence and 
the characterization of other phases occurring at higher densities. 
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Appendix A: Definition of fused hard sphere helices 

A helix is made of A^ spheres, whose centers are located at the points defined by the para- 
metric equations: 

Xi = r cos(27rtj) 

yi=r sin(27rt,) 1 < i < A^ (Al) 

^ Zi = pti 
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where r is the radius and p is the pitch of the hehx (see Fig. [I]). Given the values of r, p 
and of the contour length L, the increment At = tj+i — tj is determined by the equation: 



i „_.. ._, . /p\^ 



- = 2.AtX^+(iL). (A2) 



The Euclidean length of the helix is defined as A = zn — zi, depends on the pitch and radius, 
and coincides with the contour length L only for r = 0. 

Appendix B: Molecular volume and effective volume of helices 

The volume vo of a linear chain formed by m fused hard spheres (LHSC) of diameter D and 



center-to-center distance dec (see Fig. 12 ) is given by 

TT ^o m — 1 / „ d, 

^^0 = TT^ 



6 




(Bl) 



The same expression holds for a helix of fused hard spheres, provided that there are only 
two-sphere overlaps and the correct value of the distance d^c is used.^ For a given length of 
the curve connecting the centers of a pair of subsequent spheres, this distance depends on 
the helix radius and pitch. Table [l] reports the (geometric) volume calculated for the all the 
helices shown in Fig. |3] 



helix dec vq 

p= 2, r= 0.2 0.687 6.89 
p= 2, r= 0.4 0.680 6.85 
p= 4, r= 0.2 0.711 7.02 
p= 4, r= 0.4 0.707 7.00 
p= 8, r= 0.2 0.714 7.04 
p= 8, r= 0.4 0.714 7.04 

LHSC 0.714 7.04 



TABLE I: Volume vq of fused hard sphere helices of radius r and pitch p, calculated using Eq. 



(Bl). For comparison, also the value for LHSCs is reported. 



^ dec = \/{xi+i — Xi)-^ + (ni+i — yip + {zi+i — Zi)"^ is the distance between centres of two subsequent 
spheres. 
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A definition of the effective volume has been proposed for LHSCs, as the volume enclosed 
by the surface drawn by a sphere identical to those of the chain, rolling over the particle 



j. An example of this surface is shown in Fig 12 The effective volume of the LHSC is 
then given by the expression: 



„LHSC 



TT 



-D' 



1 + [m 




1 fd. 
2 



cc 

I) 



\ 



1 



dcc\ \ . fdcc 
2D h'^^^n^ 



(B2) 




FIG. 12: Surface defining the effective volume of a pair of fused hard spheres. 



We have adopted the same definition of the effective volume for fused hard sphere helices. 
However in this case, depending on the helix curvature, the effect of the rolling sphere can 
go beyond that of simply filling the voids between subsequent beads. We have calculated the 
effective volume (f ef ) of helices using the program MSMS [43j . The rolling sphere radius was 
taken equal to the radius of the fused hard spheres that form the helix. Table [TT] reports the 
fef values obtained for the helices shown in Fig. |3] for comparison we report in the table also 



the volume calculated according to eq. B2, using for each sphere the appropriate d^c value 
j^^LHSC-j^ ^^Q ^^^ observe that v^i = fg/^^*" for all helices with longer pitch; only for p = 2 
there is some difference, more pronounced in the case with r = 0.4. This discrepancy can be 
understood considering that these helices have grooves narrower than the sphere diameter 



D (see Fig, 12) 
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helix fcf f^f^^*" 

p= 2, r= 0.2 7.24 7.20 

p= 2, r= 0.4 7.78 7.15 

p= 4, r= 0.2 7.37 7.37 

p= 4, r= 0.4 7.34 7.34 

p= 8, r= 0.2 7.39 7.39 

p= 8, r= 0.4 7.39 7.39 

LHSC 7.39 7.39 

TABLE II: Effective volume of fused hard sphere helices of radius r and pitch p, calculated either 
by the program MSMS [43j or using Eq. (B2) with the dec distances reported in table ^ For 



comparison, also the value for the LHSC is reported. 
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